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Abstract 

Radiation-hard ionization chambers were tested using an intense electron beam from 
the accelerator test facility (ATF) at the Brookhaven National Laboratory (BNL). 
The detectors were designed to be used as the basic element for monitoring muons 
in the Main Injector Neutrino beamline (NuMI) at the Fermi National Accelerator 
Laboratory (FNAL). Measurements of linearity of response, voltage dependence, 
and the onset of ionization saturation as a function of gap voltage were performed. 



1 Introduction 



Future neutrino experiments will 
drive neutrino beams to intensity 
levels where beam secondary and 
tertiary monitoring systems must be 
concerned about effects of radiation 
damage on materials and detector 
response linearity. We have designed 
a radiation-hard ceramic pad ion- 
ization chamber (PIC) with small 
(adjustable 1-5 mm) gap thickness 
suitable for this environment. [1] 
The device will be used to moni- 
tor secondary hadrons and muons 
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in the NuMI neutrino beam for 
the MINOS experiment. This paper 
presents measurements performed 
at the ATF which demonstrate that 
the system can be used for monitor- 
ing charged particle intensities in the 
range 10 6 — 10 8 particles/cm 2 in an 
8/isec beam spill, suitable for which 
are the expected intensities at the 
MINOS muon monitors. 



2 Minos Beam Environment 



The NuMI beam, composed pri- 
marily of Up, will be produced from 
decays of secondaries generated 
when protons from the Main Injec- 
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tor at FNAL strike a graphite tar- 
get. [1] The neutrinos emerge from 
7r + — ► /U + z/^ and K + — > /i + z^ decays 
in a 675 m decay pipe beginning 
50 m downstream of a double horn 
focusing system. 

In order to predict the beam flux at 
the far detector MINOS plans to use 
the neutrino beam spectrum mea- 
sured at the near detector. However, 
neutrino events are collected too 
slowly at the near detector for spill 
by spill monitoring of the beamline 
performance. Beam monitors will be 
used to measure the distributions of 
secondary particles at the end of the 
decay pipe and of muons produced 
in the K and tt decays at three lo- 
cations (muon alcoves) in the earth, 
see Fig. 1. This system is designed to 
be sensitive to beamline component 
alignment and performance and to 
provide feedback on the beam on a 
short time scale. [1] 

In order to monitor the muons in the 
beam, arrays of PIC detectors will be 
located immediately downstream of 
the hadron absorber and successively 
after 12 m and 18 m of earth in the 
muon alcoves. To ensure that these 
chambers will operate in the linear 
region, we have measured the satu- 
ration point (defined later in Sec. 7), 
the ionization response yield, the sta- 
bility of the yield and its dependence 
on the applied electric field, as de- 
scribed in this paper. 




— 50m 675m 



Fig. 1. Schematic of the NuMI neutrino 
beamline. The absorber region and the 
locations of the beam monitoring sta- 
tions are indicated in the right hand side 
of the figure. Figure is not to scale. 

3 ATF Beam 

The Accelerator Test Facility (ATF) 
at Brookhaven National Laboratory 
provided the test beam to study the 
PIC response. The facility delivered 
42 MeV electron bunches with pulse 
lengths of ~10 ps and a range of in- 
tensities from 1 pC to 1 nC. [2] The 
beam spot size was typically between 
1.0 and 1.2 cm 2 (85% containment), 
thus the typical intensity range avail- 
able at the facility was 6x1 6 — 6x1 9 
particles/cm 2 / 10 ps. 



4 PIC Description 

The PIC is made of two metalized 
ceramic pads separated by laser cut 
ceramic washers at the corners to 
maintain the gap spacing. The ce- 
ramic pads are approximately 7.4 cm 
x 7.4 cm, a signal pad is shown in 
Fig. 2. A 0.304 mm gap separates 
the 1 cm thick guard ring from the 
active pad. The metalized-ceramic 
pad is produced commercially using 
a Platinum-Silver alloy and 1 mm 
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Fig. 2. Front view of PIC signal pad. 

thick ceramic wafers as the substrate 
by Amitron, Inc. of N. Andover, MA, 
USA. The components of the basic 
PIC unit have been chosen for their 
radiation hardness and flatness of 
less than 1.35 /iiii per cm. The radia- 
tion dose in the muon alcoves will be 
approximately 10 MRads/year. [?] 
Mechanical Properties of ceramic 
materials have been tested up to 
total radiation doses of 10 GRad. [3] 



5 Experimental Setup 

The experimental setup is shown in 
Fig. 3a. The PIC detector pads were 
mounted in a thin window gas box 
upstream of the Faraday cup (FC) 
and the beam strip chamber (BSC). 
The FC measured the total beam 
charge per pulse. The BSC was used 
to measure the beam profile and con- 
sisted of a pair of strip ionization 
chambers. Each chamber had a set of 
seven 5 mm wide strips: one set ori- 
ented horizontally and the other ver- 
tically. Each had a gap size of 2.5 mm 
per plane. Both the BSC and the 
FC completely covered the expected 
beam size of < 1.5 cm 2 (see Fig. 3b). 



Fig. 3. Schematic of the apparatus 
showing (a) the detector locations and 
(b) the front view of the beam strip 
chamber shadowed by the Faraday cup 
showing the active area of each. 

The BSC and the PIC detectors 
were operated with helium flowing 
at a rate of 100 ± 20 cc/min and at 
a pressure of 0.15 ± 0.05 inches of 
water above ambient pressure. The 
FC was also operated in a contin- 
uously purged atmosphere of dry 
helium. We used laboratory grade 
(UHP) helium independently rated 
to be 99.999% pure by Praxair Inc., 
Seattle, WA, USA. The helium gas 
was provided separately to the PIC, 
BSC and the plastic enclosure sealed 
with kapton tape around the entire 
assembly. A gas distribution mani- 
fold allowed us to control the flow 
rate and pressure to each detector. 
We used 0.25 in polyethylene tub- 
ing with brass fittings to distribute 
the gas. No provision was made to 
sample the gas for purity. 

Voltage to the detectors was pro- 
vided by two independent high volt- 
age units with local filters to suppress 
ripples. For all of the data, the BSC 
voltage remained fixed at 250 V. 
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5. 1 Electronics 



The essential features of the elec- 
tronic system used to acquire data 
are shown in the simplified schematic 
of Fig 4. An operational amplifier re- 
moved excess charge from either an 
ionization chamber with capacitance 
18 pF or the Faraday cup and the 
connecting coax cable, and stored 
it on a precision 100 pF integration 
capacitor. The voltage on the capac- 
itor was subsequently digitized by an 
ADC for computer read out. Volt- 
age on the input to the integrator 
was maintained at a virtual ground 
subject to the limitations of the am- 
plifier. This differs from pulse height 
read out systems that employ a load 
resistor to develop a voltage at this 
point in the circuit. 

The operational amplifier is part of 
a Burr-Brown ACF2101 integrated 
circuit. [4] The complete integrated 
circuit contains the necessary FET 
switches for clearing the integration 
capacitor and gating input/output 
signals. Notable limitations of the 
amplifier are a 3 V//isec slew rate 
and a random 1 mV charge transfer 
noise at the output. 

Charge integration electronics dis- 
tinguish these ionization chamber 
measurements from previous acceler- 
ator measurements which used pulse 
height read out. [5,6] Pulse height 
is more susceptible to saturation ef- 
fects because pulse height and shape 
depend on charge collection time as 
well as the amount of charge pro- 
duced in the detector. 




Fig. 4. Schematic of the ionization 
chamber and input stage of the elec- 
tronics. 

The complete system used for data 
acquisition was developed by FNAL 
for read out of a segmented wire ion- 
ization chamber (SWIC). [7] The sys- 
tem used here differs from the FNAL 
system by the addition of a bus ter- 
mination resistor to reduce channel 
cross talk for large output signals. 

A short integration gate time of 
5 msec for readout was necessary to 
prevent a small charge loss that was 
observed for much longer gate times. 

In the SWIC electronics, a large con- 
tribution to the error comes from 
random charge transfer errors of 
0.1 pC. This contributes significantly 
to the measurement uncertainty at 
FC charges of m 1 pC. All pedestals 
were taken while the beam was dis- 
abled. The pedestals did not vary 
with the PIC operating voltage. The 
pedestals for the PIC and FC show 
stability over time and show no 
signs of systematic deviation which 
would be manifest if, for example, a 
dark current were present. For the 
PIC detector, a typical pedestal was 
measured to be 0.39 pC with a typ- 
ical width of 0.13 pC. The pedestal 
width includes both statistical and 
the dominant charge transfer error. 
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6 Analysis 

We recorded the integrated charge 
for the PIC, all channels of the BSC 
and the FC for each ATF pulse. Data 
were pedestal subtracted and then 
selected for beam containment. The 
ratio of the the PIC to FC charge in 
the linear region was examined and 
anomalous events were removed. 

The mean horizontal beam position 
was determined from the BSC using 
X c = Yiqii/Yiqi where i is the strip 
number and qi is the charge collected 
on strip %. Y c , the vertical position 
was determined similarly. The hori- 
zontal and vertical beam widths were 
estimated using the standard devi- 
ation (RMS) of the beam intensity 
in each view. For the horizontal, this 



width is ax c = 




is similarly defined for the vertical 
width. 

The beam position was stable for all 
charges and the beam width was usu- 
ally less than one strip. Cuts were 
made on the beam position and hor- 
izontal and vertical widths to select 
well behaved beam pulses. During a 
typical run, the efficiency of the cuts 
was 99%. 



7 Results 

In gas ionization detectors, particles 
ionize the atoms of the gas generat- 
ing electrons and ions which in turn 
generate a detectable current if an ex- 
ternal electric field is applied. With 



a suitable electric field, the detectors 
operate in ionization mode where the 
number of ionized gas particles is sim- 
ply proportional to the total num- 
ber of incident charged particles. If 
the applied field is too large, the ion- 
ized charges are accelerated enough 
to produce additional ionization, and 
gas amplification occurs. 

Saturation occurs when the applied 
field is not large enough to over- 
come the field produced by charges 
from ionization. In this case, ionized 
charges can recombine and any fur- 
ther increase in beam intensity does 
not yield a corresponding increase 
in ionization. At voltages below the 
gas amplification region, saturation 
can be overcome by increasing the 
external electric field or decreasing 
the amount of material in the active 
ionization region so that less gas is 
ionized. The literature on saturation 
does not provide a complete physi- 
cal description under all experimen- 
tal conditions. [5,6] The saturation 
point can be experimentally deter- 
mined by measuring the deviation of 
the response from linear or by mea- 
suring the ionization response yield, 
i.e. plateau curves as a function of 
voltage. 



7. 1 Linearity 

The linearity measurements for the 
5 mm PIC at 250, 500 and 750 V are 
shown in Figs. 5, 6 and 7. The linear- 
ity measurements for the 3 mm PIC 
were performed at 250 and 400 V. 
The 3 mm 250 V data are shown in 
Fig 8; the 400 V data are similar. The 
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Fig. 5. Response of the PIC versus the 
FC charge for 5 mm 250 V data. 

results of the fits and fit ranges are 
summarized in Table 1. All data are 
shown with error bars giving the sta- 
tistical uncertainty in each bin and 
the 0.1 pC SWIC charge transfer er- 
ror added in quadrature. 

The fitting procedure [8] used the 
values of the FC and the PIC charge 
with statistical and systematic er- 
rors. For a linear fit, the data were 
binned in FC and PIC charge. 

An important input to the fit is the 
choice of fit cutoff. Since the data 
are expected to deviate from a linear 
model at charges beyond the satura- 
tion point, the x 2 /DOF was U sed to 
determine the optimal fit cutoff. We 
minimized the x 2 /DOF 

function 

of the fit cutoff. The fit cutoff, chosen 
in this way, is one objective measure 
of the deviation from linear behavior. 

All data in Fig. 5-8 are clearly lin- 
ear below the fit cutoff points, and 
saturation starts above this cut- 
off. All 5 mm ionization yields (ion 
pairs/cm) in Table 1 are consistent 
at the 1.1 a level. The fit y-intercept 
is consistent with zero. The weighted 
average ratio of the 3 mm and 5 mm 
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Fig. 6. Response of the PIC versus the 
FC Charge for 5 mm 500 V data. 
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7. Response of the PIC versus the 
Charge for 5 mm 750 V data. 



3mm BC 250V data I 




Faraday Cup Charge (pC) 

Fig. 8. Response of the PIC versus the 
FC Charge for 3 mm 250 V data. 

ionization yields is 0.60±0.01, in ex- 
cellent agreement with the expecta- 
tion of 3/5. 



6 



Table 1 

Fit parameters from linear fits to 5 mm and 3 mm gap helium data. 



Configuration 


Slope 
(pC/pC) 


PIC inter- 
cept (pC) 


X 2 /DOF 


Range 
(PC) 


5 mm He-250V 


9.95±0.22 


-0.80±1.15 


17.3/18 


0-8 


5 mm He-500V 


9.77±0.09 


0.22±0.76 


30.8/32 


0-13 


5 mm He-750V 


9.57±0.26 


1.39±1.62 


15.1/15 


0-11 


3 mm He-250V 


5.98±0.03 


1.51±0.32 


58.6/52 


0-18 


3 mm He-400V 


5.82±0.03 


2.05±0.30 


48.4/42 


0-19 



7.2 Voltage Dependence 



A voltage plateau demonstrates the 
stable operating range of an ioniza- 
tion chamber. For this measurement, 
the ATF beam was held at approx- 
imately constant intensity while the 
ionization chamber voltage was var- 
ied from a few volts up to 1000 V. 

Using the selection criteria described 
above, the average PIC to FC yield 
and RMS were determined from the 
data. The beam current at lower in- 
tensities was stable, but at higher in- 
tensities varied by ^20%. 

Plateau curves for the 5 mm and 
3 mm gap PIC in helium gas are 
shown in Fig. 9. The four curves 
correspond to beam intensities of 
5.7±0.4, 11.9±1.7, 26.8±2.4 pC for 
the 5 mm and 30.2±2.2 pC for the 
3 mm PIC. At the lowest intensity 
the ratio PIC/FC rises very quickly 
with voltage and remains constant 
beyond ~ 100 V. At the highest 
intensity PIC/FC continues to rise 
even at 1000 V indicating saturation. 
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Fig. 9. PIC response as a function of the 
gap high voltage for different nominal 
intensities for 5 mm and 3mm gap PIC. 
The lines are not a fit to the data points. 

7.3 Saturation 



To quantify saturation we plot the to- 
tal collection efficiency which is de- 
fined as the data minus the fit di- 
vided by the fit. We then character- 
ize saturation as a percent deviation 
from zero. Fig. 10 shows collection ef- 
ficiency versus FC charge measured 
at 250, 500 and 750 V for 5 mm data. 
Table 3 shows the saturation points 
as determined from the \ 2 fit (the 
fit cutoff defined earlier divided by 
the area) and 5% and 10% satura- 
tion points for the 3 mm and 5 mm 
data. We found the fit cutoff to be a 
less sensitive measure of the satura- 
tion point. The data at higher volt- 
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ages saturate at higher intensities, as 
expected. 

The ion chamber, where one pad is 
held at V(z = 0) = and the other 
pad at V(z — d) —V , is analogous 
to the parallel plate capacitor. When 
a beam of electrons passes through 
the gas, this creates a uniform charge 
distribution, p. The electric field, E, 
in this case is related to the charge 
distribution by the Poisson equation, 
dE/dz = p/e . Integrating and ap- 
plying the voltage boundary condi- 
tion, the electric field is given by, 

V ^ eo V 2 / d 

Saturation in this naive model oc- 
curs when E[z — d) — 0, which cor- 
responds to a critical charge den- 
sity p c = 2eoV/d 2 or using the ob- 
served ionization response yield of 
19.5 ion pairs/cm and d = 0.5 cm, 
p c = 0.708 x V pC/cm 3 , with V in 
volts. 

To calculate the intensity per unit 
area the area of the beam spot on the 
chamber was measured using three 
different methods. For each method, 
we used the BSC data averaged in 
the linear region, which always ex- 
tended to intensities larger (due to 
the smaller gap of the BSC) than 
the saturation point of the PIC. The 
beam spot was observed periodically 
(for every fifth pulse) with a fluo- 
rescent screen viewed with a CCD 
camera and found to be approxi- 
mately described by an ellipse. The 
area given by, 

A = 7ia Xc (T Yc , 



with major and minor axes 2ax c 
and 2<7y c was determined by fitting 
a gaussian to the BSC profile for the 
85% contained beam area. Fig. 11 
shows that the measurement of the 
the area was stable and did not de- 
grade as a function of intensity. 

Two additional methods made use of 
the segmentation of the BSC. A grid 
of values was formed by the outer 
product of the two vectors consisting 
of the charge measured in the hori- 
zontal and vertical strips. The result- 
ing matrix was then normalized by 
the total charge. The first method, 
("Box I"), simply added the number 
of grid elements with more than 5% 
of the maximum charge. The second, 
("Box II") summed each grid area 
weighted by its normalized charge. 

The results are shown in Table 2. 
The three methods agree within 20% 
for all measurements. The average of 
these three measurements is taken 
to be our area estimate. The errors 
on the measurements are highly cor- 
related and the measurements them- 
selves are consistent. Thus, the error 
on the average was chosen to be the 
largest percent error from the in- 
dividual measurements. Due to the 
large size of the BSC strips compared 
to the beam size, both the area and 
error are likely overestimated. 

Fig. 12 shows the intensity at which 
the data saturate as a function of 
V/d 2 for the 5% and 10% level. The 
linear fits shown in the figure in- 
clude only the 5 mm data. The 3 
mm data are in good agreement 
with the fit except for the 400 V 
10% data point which is lower by 
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Table 2 

Area estimates for 5 mm and 3 mm data. 



Data 


Method (cm 2 ) 




TT" /T t 7- /Tt r 

"° X C °Y C 

(85% Containment) 


R nY T 


Rnv TT 


5 mm 250 V 


1.16±0.14 


1.07±0.14 


1.15±0.10 


5 mm 500 V 


1.39±0.13 


1.55±0.21 


1.08±0.15 


5 mm 750 V 


1.16 ± 0.17 


1.26 ± 0.25 


1.12 ± 0.13 


3 mm 250 V 


0.99±0.25 


1.34±0.22 


0.91 ± 0.21 


3 mm 400 V 


0.88±0.11 


1.26±0.16 


0.81 ± 0.15 
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Fig. 10. PIC charge collection efficiency 
versus integrated beam charge mea- 
sured in the Faraday cup for the 5 mm 
data. 
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Fig. 11. The 85% contained area plots 
show the area as a function of intensity 
for 5 mm 250 V data. 

2.5 sigma. If the 3 mm data are in- 
cluded in the fits, the fit to the 5 
% data points does not change sub- 



stantially, but the x 2 /DOF for the 
10 % fit worsens from 1.1 to 2.1. The 
5 mm data demonstrate the naive 
model expectation that the satura- 
tion point behaves linearly in V/d 2 . 
For the 5 mm 500 V data, the crit- 
ical charge density corresponds to 
11.3 x 10 7 particles/cm 2 , which is 
in good agreement with the mea- 
sured value of saturation 11± 2 xlO 7 
particles/cm 2 for the 5% saturation 
point. We note that we have not at- 
tempted to understand the effects 
of multiplication on ionization yield 
as well as on the apparent satura- 
tion point. At high electric fields and 
intensities effects of multiplication 
could make the interpretation of the 
data more complex. 



8 Conclusions 

Results from the ATF measurements 
have demonstrated that the PIC de- 
vice with a 5 mm or smaller gap 
spacing is suitable for operation in 
the high-rate environment expected 
at the muon monitors of the NuMI 
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Table 3 

Saturation limits for 5 mm and 3 mm data. 





Saturation Measure 




X 2 Fit 


5% Point 


10% Point 




10 7 particles/cm 2 


o mm zoU V Data 


4.4 ± 0.6 


5.0 ± 0.7 


9.4 ± 1.2 


5 mm 500 V Data 


6.0 ± 0.8 


11 ± 2 


22 ± 3 


5 mm 750 V Data 


5.8±1.0 


14 ± 3 


23± 4 


3 mm 250 V Data 


10 ±3 


18 ± 5 


23 ± 6 


3 mm 400 V Data 


12 ±3 


20± 4 


23 ± 5 
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Fig. 12. Saturation point as a function 
of the voltage. Data for the 5% and 10% 
saturation points are shown. The two 
linear fits include only the 5 mm data. 



beam. The 5 mm PIC operated lin- 
early with stable response over a 
voltage range 250-750 V in helium 
gas with intensities up to 14 (23) 
xlO 7 particles/cm 2 for a 5% (10%) 
saturation for a short-pulse beam. It 
is expected that for longer duration 
beam spills of order 10 /is, the satu- 
ration point should increase and this 
result will be a lower bound. Further- 
more, the results show that satura- 
tion can be modified in a controlled 
manner by increasing or decreasing 
the externally applied field. 
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